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Relativistic current carrying strings moving axisymmetrically on the background of a Kerr black 
hole are studied. The boundaries and possible types of motion of a string with a given energy 
and current are found. Regions of parameters for which the string falls into the black hole, or is 
trapped in a toroidal volume, or can escape to infinity, are identified, and representative trajectories 
are examined by numerical integration, illustrating various interesting behaviors. In particular, we 
find that a string can start out at rest near the equatorial plane and, after bouncing around, be 
ejected out along the axis, some of its internal (elastic or rotational kinetic) energy having been 
transformed into translational kinetic energy. The resulting velocity can be an order unity fraction 
of the speed of light. This process results from the presence of an outer tension barrier and an inner 
angular momentum barrier that are deformed by the gravitational field. We speculatively discuss 
the possible astrophysical significance of this mechanism as a means of launching a collimated jet 
of MHD plasma flux tubes along the spin axis of a gravitating system fed by an accretion disk. 
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I. INTRODUCTION 

In this paper we study the motion of an axisymmetric, 
current carrying relativistic string loop in the background 
of a Kerr black hole. This interesting dynamical system 
has been studied before in its own right [1, 2]), and our 
goal here is to better understand its behavior. 

Our investigation is loosely motivated by the prob- 
lem of the production of collimated astrophysical jets of 
matter. Systems that exhibit collimated jets range from 
accreting young stars, neutron stars and black holes to 
supermassive black holes in quasars and active galactic 
nuclei [3] . Magnetized plasmas interacting with an accre- 
tion disk are believed to play a central role in jet produc- 
tion, but despite much work and many prosposals for the 
mechanism, the process is still not yet well understood. 
Such plasmas are governed by magnetohydrodynamics 
(MHD), a complicated nonlinear field theory. Under 
some circumstances, plasmas exhibit associated string- 
like behavior, cither via the dynamics of magnetic field 
lines embedded in the plasma [4-6] , or in the dynamics of 
relatively thin isolated flux tubes of plasma which can be 
approximately described by an effective one-dimensional 
string [7]. In cither of these cases, the string would be 
described by an energy density and a tension, and might 
be carrying currents of mass and/or charge. The hope is 
that some essential aspects of the physics could be cap- 
tured by string dynamics, which is tremendously simpler 
than MHD. 

The dynamical formalism for describing relativistic 
current carrying strings is well-developed [8-11], and a 
generalization to a wider variety of equations of state is 
given in Ref. [12]. Here we restrict attention to a rela- 
tively simple class of strings, namely axisymmetric loops 
characterized by a constant tension and a conserved mass 
current on the string worldsheet, in order to begin devel- 
oping some insight into the factors influencing the string 
dynamics. This could be generalized in future work to 
allow for different equations of state, coupling of a charge 



current to an external electromagnetic field [13, 14], and 
deviations from axisymmetry [15-18]. 

The string tension prevents a string loop from expand- 
ing beyond some radius, and the worldsheet current can 
produce an angular momentum barrier preventing the 
loop from collapsing into the black hole. Thus, depend- 
ing on its energy and current, a string may be trapped in 
a toroidal region surrounding the black hole, or its mo- 
tion may be confined to a cylindrical shell that extends 
to infinity. In the latter case, one can envisage processes 
wherein internal energy of the string is converted into 
translational kinetic energy via the propagation on the 
black hole background. This is one of the phenomena we 
aim to understand. The other key question is what role 
may be played by the spin of the black hole, as a result 
of the associated dragging of inertial frames. 

II. CURRENT CARRYING STRING IN 
CURVED SPACETIME 

The string worldsheet is described by giving its space- 
time coordinates X x {a a ) (A = 0,1,2,3) as functions of 
the two worldsheet coordinates a a (a = 0,1). The in- 
duced metric on the worldsheet is 

h ab = 9»»X%X% (1) 

where g^ is the ambient spacetime metric, for which we 
adopt the signature ( — h++). To describe the current we 
introduce a scalar field (p(<r a ) living on the worldsheet. 
We consider dynamics generated by the action 

S = -Jd 2 a V^h (n/c + h ab V , aV , b ) , (2) 

where the constant fx denotes the string tension, and c is 
the speed of light which is hereafter set to unity. We as- 
sign the line element g^dx^dx" the dimensions of length 
squared, and take the worldsheet coordinates to be di- 
mensionless, so that the dimensions of h a b are length 
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squared, those of \i are energy over length, and those 
of ip are square root of action. 1 

The action is the integral of a worldsheet scalar den- 
sity, hence it is independent of the choice of worldsheet 
coordinates. Note that the part of the action involving 
the scalar field ip is invariant under conformal rcscalings 
of the metric. One could of course consider other scalar 
functions of the invariant h ab p :a p :b in the Lagrangian, 
but for the exploratory purposes of the present investi- 
gation the simple choice in (2) is sufficient. In any case, 
one ultimately wants to generalize the system to allow 
for a wider class of equations of state that are perhaps 
more suitable for MHD applications. 

Varying the action with respect to h ab yields the world- 
sheet stress-energy tensor density S a6 , 

5 h S = \ $ <P<rt ab 5h ab , (3) 

where 

t ab = V^h(2.fj b - ( M + j 2 )h ab ) , (4) 

and 

3a = <P,a, f = h ab j b , f h-'\i,,il,. (5) 

The tilde on S°* is a reminder that this quantity has 
density weight one with respect to worldsheet coordinate 
transformations. Because of the conformal invariance, 
the j-dependent part of T, ab is traceless. The contribution 
from the string tension is proportional to the metric, and 
with fj, > has a positive energy density and an opposite, 
negative pressure, i.e. a tension. The contribution from 
the current is traceless, due to the conformal invariance 
of that part of the action. It can be viewed as a 1+1 
dimensional massless radiation fluid, with positive energy 
density and equal pressure. 

Varying the action with respect to X^ yields the equa- 
tion of motion 2 

(£ ab g, x X%), b - |E«V,A*,W = 0- (6) 

Varying the action with respect to ip yields the 1+1 di- 
mensional wave equation, 

(V^hh ab p, a ), b = 0, (7) 

which also implies that the current is divergcncclcss, 
W—hj a ). a = 0. The worldsheet stress-energy tensor is 
also divergenceless, T, ab - b — 0, with respect to the h ab - 
compatible covariant derivative. (To verify this directly 
one must use j a ,b = jb,a, which follows from j a = tp a . 
Conversely, (v 7 — hj a ), a — and T, ab - b = 0, together im- 
ply that j a = ip >a for some scalar field p.) 



1 Strings described by action (2) were originally introduced [19] 
as an effective description of "superconducting strings", a type 
of topological defect that might occur in a theory with multiple 
scalar fields undergoing spontaneous symmetry breaking. 

2 One could rewrite eq. (6) in terms of the Christoffel symbols of 
9liv, however it is more convenient to work with the explicit form 
(6). 



III. AXISYMMETRIC DYNAMICS IN A 
STATIONARY, AXISYMMETRIC SPACETIME 

A. Conformal gauge 

Every two dimensional metric is conformal to a locally 
flat metric, hence in particular 

h ab = fl 2 T] ab , (8) 

where rj a b is locally flat and f2 is a worldsheet scalar func- 
tion. We can always adopt coordinates a a — (r, a) such 
that the coordinate components of r\ ab are r\ Ta = and 
i] TT = —r} aa = —1, i.e. the constant Minkowski metric. 
This "conformal gauge" choice is equivalent to requiring 
h Ta = and h TT + h aa = 0, which as conditions on 
read 

9llv X%X»„ = = g„ v (X»X» T + X%X^). (9) 

In this gauge h = — ft 4 , and the metric components 
satisfy y '—hh ab = rj ab , where r/ ab is the inverse of the 
Minkowski metric. The conformal gauge conditions (9) 
do not completely fix the coordinates: the lightlikc com- 
binations r±(T may be replaced by any smooth monotonic 
functions of themselves. 

We aim to study axisymmetric string motion in an 
axisymmetric, stationary spacetimc. Such a spacetime is 
described by a metric, written in coordinates (t,r,8,(j)), 
of the general form 

ds 2 = g tt dt 2 + 2g t(f> dtd<p + g^dcf) 2 + g rr dr 2 + g ge d6 2 (10) 

where all the metric: components are independent of t and 
tj>. 

It is tempting to choose the worldsheet coordinate a 
to be equal to the spacetime coordinate <f), in which case 
a general axisymmetric worldsheet would take the form 
X^(a, t) = (t(T),r(T),8(T),a). However, this would im- 
ply h Ta — gt(f>t^, which does not satisfy the conformal 
gauge condition if g t <f, ^ 0. We therefore allow for a r 
dependent relative shift 4> — & + f( T ), so that 

X^(a,r) = (t(r),r(r),e(T),a + f(r)). (11) 

The gauge condition h Ta = then becomes 

S^"=W + fl#/ = 0. (12) 
which determines / via 

f = -(9t<t>/9tt)i, (13) 

where the dot stands for d/dr. Equation (12) means that 
X^, and thus the worldsheet vector d T , are orthogonal to 
the axial rotation Killing vector d^. That is, the curves of 
constant a are zero angular momentum (ZAMO) world- 
lines. For later reference we note that this same gauge 
condition yields the useful relation 

g^ t X^ = {g tt - g 2 t4> lg^)i. (14) 
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The second gauge condition h TT + h aa = becomes 



(gtt - g 2 t4> /g u )i 2 + g rr r 2 + g 8e 6 2 + . 



0. 



(15) 



The conformal factor ft 2 is determined in this gauge by 



the equation h„ 



fl 2 rj aa = fl 2 , i.e. 
^ 2 = g<p<p- 

B. Equations of motion 



(16) 



In the conformal gauge (8) the equation of motion (7) 
for f takes the simple form 



<P,T 



0. 



(17) 



The general solution is ip = f+(r + a) + <p-(r — <j), 
where ip± are arbitrary functions. On the other hand, 
the assumption of axisymmetry implies that the current 
is independent of a. Thus j 0)(T = 0, which implies (p, Ta — 
tp.oo = 0. Together with (17) this implies that ip must 
be linear in both a and r, 



<P = 3<t<t + 3tT, 

where j a and j T are constants. 3 Then j 2 = 
and the components of T, ab (4) become 

= 9^(jl+3 2 r) + » 



(18) 

Sl-\jl-3 2 r), 

(19) 
(20) 
(21) 



The quantity S TT in the conformal gauge is equal to 
the string energy density as measured by observers co- 
moving with the string at constant a. These are the 
co-moving ZAMO observers mentioned in the discussion 
following cq. (13). To verify the claim about the relation 
between S TT and the string energy density, note that 
since the 2-velocity u a of those observers is orthogonal 
to the constant r surfaces, we have u a = Qd a T, the pref- 
actor being determined by the unit normalization of 
u a . Also h = —ft 4 . The energy density is thus given by 
(— /i)~ 1 / 2 S afc u a it; ) = S TT . In particular, /i is the contri- 
bution form the string tension, and (j 2 + j 2 )/g<p<p is the 
energy density of the current. 

The equations of motion (6) can now be written out 
more explicitly in our chosen gauge. Using the gauge 
conditions and the axisymmetry, they take the form 



(£ rr 5AlA A^ + £^ A ), T 



(22) 

, A + £ TCT g t 0,A i) = 0. 



3 Note that while j s is single-valued, the scalar ip is multi-valued. 
When the string action is taken to be an effective description of 
a field theory defect [19], <p is a phase defined modulo 2tt, which 
implies that the possible values of the current are quantized. 



For A equal to (p this equation is identically satisfied. For 
X = t, it yields (using (14)) the energy conservation law, 



Z TT (gtt - 9%lgu)i + ^ aT 9t4> = -E, 



(23) 



where E is a constant to be identified below with the 
total Killing energy of the string (divided by 2tt). For A 
equal to r or 9 it yields 

(t^g rr r) ^ - (t TT g^ r X»X v + t™g^. r ) = (24) 

(± TT gee0) ^ - (± TT g^ e X»X v + g^,e) - 0. (25) 

The dynamics depends on the current only through 
the worldsheet stress tensor (19-21), in which the current 
enters only in the two quadratic combinations j 2 +j 2 and 
jcrjr , both of which are symmetric under interchange of 
ju and j T . Note that j a j T appears only via t, through 
the energy conservation law (23), and then only when 
5t0 + 0. 

To parametrize the solutions, we will later use the first 
of these combinations, together with the ratio of the cur- 
rent components, 
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-ja/jr- 



(26) 

(27) 



(The minus sign is included so that positive to will 
correspond to positive angular momentum. Note that 
—jo/jr = 3° / j T ■) The product j a j T can be expressed 
in terms of J and w as j a j T = —ujJ 2 /(\ + oj 2 ). Since 
the dynamics is symmetric under interchange of j a and 
j T , the range —1 < u < 1 covers all distinct cases. The 
extreme values u> = ±1 corrcsond to null currents. 



C. Conserved quantities 

Suppose £ A is a spacctimc Killing vector. If the space- 
time coordinates are chosen so the components of £ A 
are constant everywhere, then Killing's equation implies 
gnu,\£, X — 0. Contracting (6) with £ A then yields a con- 
served Killing current. That is, the worldsheet vector 
density 



31 = ^ ab x% 9flX e 



satisfies 



j£ ,b 



0. 



(28) 



(29) 



Although we argued for its conservation using an adapted 
coordinate system, is manifestly spacetime coordinate 
independent. Note that it is just the worldsheet energy- 
momentum tensor contracted with the pullback of the 
Killing one-form to the worldsheet. Thus integrating it 
over any closed cross section of the worldsheet gives the 
conserved quantity associated with the Killing vector, 



(30) 



Qi = / J£dS b . 
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If we take the integral over a surface of constant r and 
use the axisymmctry, the integral becomes simply 



(31) 



When £ A is the Killing vector <9 t corresponding to t 
translation symmetry then (11) and (14) yield 



-Q dt = £ = 2ttE, 



(32) 



where E is the constant introduced in cq. (23). Thus 
E is in fact the Killing energy of the string divided by 
2ir . When £ A is the Killing vector corresponding to 
(f> translation symmetry (rotation), then using (11), (12) 
and (21) we find the angular momentum L of the string, 



-4irj a j T . 



(33) 



This is manifestly constant for the solutions (18) under 
consideration. Without the current, the string carries 
no angular momentum, since the stress tensor is then 
proportional to h a b, which is boost invariant along the 
worldsheet. 



D. Effective potential 

We can solve (23) for i and substitute into the gauge 
condition (15), yielding 

(g^/g^ - gtt)& TT )\grrr 2 + g ee 2 ) + V(r, 6) = 0, (34) 
where 

V(r, 6) = -(E + g t4> t° T ) 2 + {g 2 t4> - g tt g^)(± rT ) 2 (35) 

is what we call here the "effective potential". The mo- 
tion outside the horizon is confined to the region where 
V(r, 6) < 0, since the first term in (34) is non-negative 
outside the horizon. 4 More explicitly, the motion is 
bounded by the contour 5 



E = E b = (^-M#) 1/2 S TT -flt^ ffT (36) 
= A{r,0)p + B{r,6,Lo)J 2 , (37) 



where A and B are functions determined by the metric 
components and, in the case of B, the value of u. The 
initial conditions and the value of the current determine 
E through (15) and (23), and the subsequent motion is 
confined to the region contained within the contour (36). 
We shall determine that motion in several representative 
cases by numerical integration. 



IV. NEWTONIAN ELASTIC RING 

Before presenting the results for the relativistic string, 
we develop a Newtonian analogue of the system, to help 
interpret the dynamics. Indeed, it turns out that most of 
the salient features are independent of relativistic effects. 

Consider a Newtonian elastic ring of total mass to, 
moving axisymmetrically in the gravitational field of a 
point mass M. Using cylindrical coordinates (p,z,<f>), 
the conserved Newtonian energy of the system is given 

by 



E 



N 



\m{p 2 + z 2 + p 2 uj 2 ) + \kp 2 



GMm/(p 2 + z 2 ) 1 / 2 , 
(38) 

where u> is the angular velocity and k the elastic constant. 
The angular momentum I = mp 2 uo about the axis is con- 
served. Using I, energy conservation can be expressed 
as 



\m{p 2 + z 2 ) + V N {p,z) = 0, 



with 

V N (p, z) 



-En - 



2mp 2 



\kp 2 



GMr 



(39) 



(40) 



The motion of the ring is confined to the region 
Vn(p,z) < 0, and thus bounded by the contour 
Vn(p, z) = 0, i.e. the curve in the p-z plane given by 



E 



I 2 



N 



2mp- 



+ \kp 2 



GMm 



Vp 2 



+ z 2 



(41) 



The bounding contour is given in Fig. 1 for different 
values of the angular momentum and energy, in units 
with to = k = GM = 1. The angular momentum term 
cx p~ 2 dominates at small p and makes an inner p wall, 
while the elastic energy term cx p 2 dominates at large p 
and makes an outer p. The gravitational potential well 
deforms the inner wall inward toward the center, and the 
outer wall outward. The figure illustrates that, depend- 
ing on the energy and angular momentum, the ring can 
either oscillate in a confined toroidal region around the 
equatorial plane, or escape to infinity confined to a chan- 
nel along the azimuthal direction. In the exceptional case 
of zero angular momentum, the ring can be confined to 
a spherical region including the origin. 



V. STRING DYNAMICS 
A. String in flat background 



4 The determinant of the metric (10) is (gttg^tfy — g^^grrgee, which 
must be negative for a Lorentzian metric. For the metrics we con- 
sider, gee > everywhere and g r r,g<t>^ > outside the horizon, 
so evidently gf^/g^^ — gtt is positive outside the horizon. 

5 To obtain this we have taken a square root, but the root is unique 
since, according to (23) the combination E + gt^S CTT is positive 
for trajectories outside the horizon. 



In order to demonstrate the effect of the current on 
the motion of the string in a simplified setting, and to 
obtain a useful formula for the string energy far from 
the black hole, let us first consider the case in which 
the background is just flat spacetime. We use cylindrical 
coordinates (t, p, z, <f>), so g^ v — diag(— 1, 1, 1, p 2 ), and we 
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Figure 1: Newtonian elastic spinning ring of total mass 
m, moving axisymmetrically in the gravitational field of a 
point mass M. Contour plots of the boundaries on the mo- 
tion for different energies and for angular momentum I = 
0, 0.12, 1.5, 2 (in units with m = k = GM = 1). Both bound 
trajectories and escape trajectories can be present, depending 
on the energy. 



consider axisymmetric motion. The effective potential 
(35) then takes the form 



V(p) = -E 2 + p 2 (Z Tr ) 2 . 



(42) 



The motion is confined to the region where V(p) < 0. 
The bounds of the motion are denned by V(p) = 0, i.e. 



E = E b = pt TT = np+ J 2 /p. 



(43) 



The right hand side is l/2n times the Killing energy of 
the string when r = 9 = 0. Since the string is then at 
rest, and there is no redshift factor in flat spacetime, this 
is the same as the energy measured in the rest frame of 
the string. The first term represents the elastic energy 
of the string, while the second term is the angular mo- 
mentum barrier arising from the current circulating in 
both directions around the string. Alternatively, it can 
be viewed as a variable "rest mass" of the string, aris- 
ing from the energy density J 2 /p 2 of the current. In the 
presence of the current, the string oscillates between two 
p values, the roots of (43), never collapsing to a point. 

The elastic energy and angular momentum play similar 
roles in the relativistic case (43) and the Newtonian case 
(41), but with different powers of p. The relativistic elas- 
tic energy is proportional to the string length rather than 
its square, and the angular momentum barrier comes 



from a kinetic energy that is linear rather than quadratic 
in the momentum. A more important difference between 
the Newtonian elastic ring and the relativistic string is 
that in the latter the current can circulate in cither di- 
rection. Both directions can simultaneously contribute 
to the angular momentum barrier, so the barrier is not 
determined by the net angular momentum. Indeed it is 
present even when the net angular momentum vanishes, 
which is the case when either j a or j T vanishes. 



B. String in Kerr background 

1. General setting 

Next we consider the full problem, i.e. the motion of 
the string in Kerr spacetime, with line element (in Boyer- 
Lindquist coordinates) 



ds 1 



2Mr\ 2 4Mrasin 2 



(It'- 



ll 2 



-dtdcj) 



,2,2, 2Mra • 2q\ ■ 2q,, 2 

+ | r 2 + a A H ^ — sm I sm 



+ ^dr 2 + R 2 d9 2 . 



(44) 



Here M and a = J/M are the mass and the specific 
angular momentum of the black hole respectively, 

R 2 = r 2 + a 2 cos 2 9, A 2 = r 2 - 2Mr + a 2 , (45) 

and we are using units in which G = c = 1. (Below 
we shall instead choose units with GM = 1.) The event 
horizon is located where A = 0, at 



r h = M + \/M 2 -a 2 , 



(46) 



and the boundary of the ergoregion lies where the coeffi- 
cient of dt 2 vanishes, at 



r e = M+ ^JM 2 - a 2 cos 2 9. 



(47) 



We consider only motion outside the horizon, so this co- 
ordinate system is adequate. 

The explicit form of the effective potential (35) is 
too complicated to be illuminating, unlike the flat back- 
ground case. A given orbit is determined once j s , jt, 
and initial conditions are specified. The initial condi- 
tions determine the values of the energy E and angular 
momentum L, and the subsequent motion is confined to 
the region contained in the contour (36). The equations 
of motion can be integrated numerically to explore the 
detailed behavior of the string motion. 

We will present here only results for the string mo- 
tion around a rapidly spinning black hole, with a — 0.99. 
This is complementary to the a = case presented by 
Larsen [1]. The role of the black hole spin is addressed 
below in subsection VB5. We consider just two sets of 
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Figure 2: Graphical representation of the initial conditions, in 
units with GM/c 2 — 1. The area inside the horizon is covered 
in grey, the dotted red line marks the equatorial plane and 
the dashed black curve is the boundary of the ergoregion. Pi 
corresponds to r = ru + 1-25 and P2 to r = nsco- At both 
points, 8 — 7r/2 + 0.02, and we choose r — 6 = 0. The specific 
angular momentum is a — 0.99 and, therefore, P2 lies within 
the ergoregion (nsco <r e ). 

initial conditions for several different values of J 2 and w, 
which appear to be adequate for understanding the fea- 
tures of the motion. For both sets the string starts out 
at rest, ro = #0 = 0. The initial polar angle is fixed in 
both cases at #0 = 7r /2 + 0.02, just below the equatorial 
plane. The two cases differ only in the initial value of r. 
The two values we will explore are (somewhat arbitrar- 
ily) r a = r h +1.25M « 2.39 M and r = n SC o ~ 1.45 M, 
where nsco i s the radius of the innermost stable co- 
rotating circular orbit (ISCO) which is given in terms 
of M and a in, e.g., Ref. [20]. 

In Fig. 2 the two sets initial conditions correspond to 
points P\ and Pi respectively. Each point in this figure 
other than those on the y axis corresponds to a circular 
ring in space. The coordinates x and y, used as labels in 
this and later graphs of the (r, 8) plane, are defined by 

x = rsin8, (48) 
y = r cos 8. (49) 

The grey area indicates the region inside the horizon, 
the dotted red line marks the equatorial plane and the 
dashed black curve is the boundary of the ergoregion. 
For lo we consider only the three values lo — 0, ±1, since 
these span the full dynamical range and can be expected 
to frame the full picture. Note that changing the sign of 
lo is equivalent to changing the sign of a, since lo appears 
always multiplied by a (via g tl p) in the equations of mo- 
tion, and all other occurrences of a are via a 2 . Positive lo 



corresponds to co-rotation, while negative lo corresponds 
to counter-rotation. 



2. Boundaries of motion 

By studying the boundaries of the motion we can clas- 
sify the types of dynamics and determine how they de- 
pend on the current, energy, and initial conditions. These 
boundaries are determined by eq. (36). For a given ini- 
tial condition and current, a particular boundary is de- 
termined. These boundaries are the relativistic version 
of the contours plotted in Fig. 1 for a Newtonian elastic 
spinning ring and, as we will see shortly, the share many 
qualitative characteristics with their Newtonian counter- 
parts. At the end of this subsection we give examples of 
such boundaries, but first we present some graphs, de- 
duced from the potential, that allow one to see which 
ranges of current or energy will correspond to different 
qualitative types of boundaries. 

To begin with, in Fig. 3 we have plotted Eb as a func- 
tion of J = y/ H + j 2 and r for fixed = it/2 + 0.02, 
which corresponds to our initial conditions, and for three 
different values of lo. In this and all subsequent plots we 
adopt units with GM = c = fi = 1. That is, in terms of 
the gravitational radius r g = GM/c 2 , what is plotted is 
Eb//J*r g and J yj cjjx/r g vs. r/r g . 

In the cases lo = 0,-1, where there is a maximum out- 
side the horizon for any given value of J, the motion has 
an outer boundary in r but no inner boundary if either 
E < E™ in or E > £^ ax , where E™ in and £™ ax corre- 
spond to the minimum and the maximum of Eb respec- 
tively. If E™ ax > E > E™ in then there will be three solu- 
tions to eq. (36), corresponding to an inner and an outer 
boundary and yet another outer boundary at low values 
of r. In the co-rotating case lo — 1, where there is no 
maximum, the role of £;™ ax is played by E b , the value of 
Ef, at r = Th- For both an inner and an outer boundary to 
exist, the energy must fall in the range E' b l > E > E b nin . 
Notice that in this case no string can exist outside the 
horizon with energy E < E b mn . 

So far we have been discussing bounds in the radial di- 
rection near the equatorial plane. There is also a bound- 
ary in the motion along the axis if the total energy of 
the string is less than the energy of a string at rest at 
infinity, since the string must have at least this much en- 
ergy in order to be able to escape. We call this energy 
E CBC . To determine how it depends on J, we can use the 
flat background example of section VA. Eqn. (43) gives 
Eb = [ip+J 2 /p, whose minimum lies at E CSC (J) — 2Jy/J2. 

Fig. 4 presents plots that allow one to see how the 
qualitative properties of the motion depend on J, for 
specific values of lo and the initial r as indicated, and 
for 8 = it/2 + 0.02. The solid black curve, which can be 
thought of as constant r slice of the graphs in Fig. 3, is 
the energy (36) of the string at rest at the given radius 
r, and is the only curve on the graph which depends on 
the value of the initial r. The dependence on lo comes in 



7 




(a) u) = 1 




(b) a; = 
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(c) u = -1 

Figure 3: Plots of the energy of a string at rest (or energy at 
the boundary of motion) as a function of r (plots starting at 
r = Vh) and + for three values of the ratio between 
the components of the current u> = —ja/jr and at 9 — n/2 + 
0.02 (GM = c = n = 1). All graphs possess minima, which 
implies that bound motion is possible for energies greater than 
-E™ ln , the energy associated with the minimum. However, 
only for u) = 0,-1 is there a maximum located outside the 
horizon. 



only through the g tc j, term in (36). Increasing to acts to 
raise the energy for a given ,/. The dashed red line cor- 
responds to the energy of a string at rest at infinity (the 
escape energy). The shaded region is the area between 
the curves corresponding to the minimum and the max- 
imum of Ej, for each J (cf. Fig. 3), with the exception 
of the ijj = 1 cases where the dotted-dashed black curve 



corresponds to the value of £), at the horizon. 6 

We now demonstrate how one can infer all the qual- 
itative information about the boundaries of the motion 
from Fig. 4. For a given energy or value of J, there is a 
corresponding point on the solid black curve. This point 
can lie in several distinct regions of the graph, which cor- 
respond to qualitatively different types of boundary: 

(i) A region outside the shaded area and above the 
dashed red line. In this case there is just an outer 
boundary, i.e. there is only one value of r outside 
the horizon that yields a solution to eq. (36) for each 
for given E and J. The string can potentially es- 
cape in the y direction, as its energy is higher than 
the E csc . 

(ii) A region outside the shaded area and below the 
dashed red line. Again, there is just an outer bound- 
ary. However, now the string cannot escape in the 
y direction as its energy is lower than the -E e sc- 

(iii) A region inside the shaded area and above the 
dashed red line. In this case there are always outer 
and inner boundaries, and for the examples with 
u> = 0, — 1 there is yet another outer boundary at a 
lower value of x. That is, eq. (36) has up to three 
solutions for r outside the horizon each value of 6 
and for given E and J. However, there is no bound- 
ary on the y motion and the string can potentially 
escape. 

(iv) A region inside the shaded area and below the 
dashed red line. The boundaries of the motion 
are the same in number as in the previous case, 
but with a qualitative difference: the inner and the 
outer boundaries at higher values of x meet at a 
certain value of y and together form a closed curve. 



It is rather surprising to have a string (instantaneously) at rest 
at the horizon with non-zero Killing energy. This comes about 
as a singular limit whose nature is hidden by a singular gauge. 
The first term in eq. (36) vanishes on the horizon, so we have 
there E b = -g t<#) E <TT = 2(g t ^/ g^jajr , using eq. (21) in the 
last step. The metric factors are finite and non-zero as long 
as the black hole is spinning, in which case the energy is non- 
zero as long as both j a and j T are non-zero. But if the string 
is instantaneously static at the horizon, then the worldsheet is 
null on that slice, rather than being time-like, so the conformal 
gauge is not accessible. In the limit as such a configuration is 
approached, the conformal gauge tangent vector X^ T approaches 
the null horizon generator, becoming infinitely stretched since 
the gauge condition requires that it has a nonzero timelike norm. 
Thus the current 4- vector jt 1 = X%j T + X^j" itself diverges as 
long as j T = h TT j T = g^jr is non-zero. A non-zero Killing 
energy configuration of an instantaneously static string at the 
horizon arises from this singular limit of the current, but only, 
as seen above, provided the black hole has spin. The underlying 
reason for this last requirement is that only in the presence of 
spin is the Killing vector dt, with respect to which the energy is 
defined, distinct from the null Killing vector that is normal to 
the horizon. 
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Figure 4: Plots of the energy of a string at rest (or at the boundary of motion) as a function of J = \J j% + j| for specific values 
of ui and r as indicated, and for 6 — 7r/2 + 0.02 [GM = c = fj, = 1). The solid black curve, which can be thought of as a constant 
r slice of the graphs in Fig. 3, is the energy of the string at rest at the given radius r, and is the only curve on the graph that 
depends on the value of r. The dashed red line corresponds to the energy of a string at rest at infinity (the escape energy). The 
shaded region is the area between the curves corresponding to the minimum and the maximum of for each J (cf. Fig. 3), 
with the exception of the u = 1 cases where the dotted-dashed black curve corresponds to the value of Ef, at the horizon. These 
graphs allow one to infer the properties of the motion. For a given energy or value of J, there is a corresponding point on the 
solid black curve. If this point lies outside the shaded area and above the dashed red line then the motion has just an outer 
boundary in the x direction, but the string has enough energy to potentially escape (there is no boundary in the y direction). 
If the point lies outside the shaded area and below the dashed red line then the motion is bounded in all directions. If instead 
the point lies within the shaded region, then the motion has boundaries both at low and high values of x. Additionally, if the 
point lies above the dashed red line then there is no boundary on the y motion and the string can potentially escape. If the 
point lies below the dashed red line then the motion can be trapped in some toroidal volume. Notice that for ui = — 1 and 
inside the ergosphere, Fig. 4(f), there can never be an inner boundary. 



Therefore, there is also effectively a boundary in 
the y motion and the string can be trapped in some 
toroidal volume. 

There is a subtle point that should not pass unnoticed: 
depending on the initial conditions, a string with given E 
and J which corresponds to a point in the shaded region 
in Fig. 4 can either be moving in between the inner and 
outer boundary at higher it's, or it can be inside the 
second outer boundary at smaller x's. Finally it should 
be stressed that a second outer boundary at small x's 



can only exist when there Eb has a maximum outside 
the horizon. 

Next we give in Fig. 5 a few examples of plots of the 
boundaries of motion corresponding to given initial con- 
ditions and current (or energy), for r = r/j + 1.25. The 
qualitative nature of these boundary plots can be read 
off from Fig. 4. For example, with E = 50 there exist 
both inner and outer boundaries in the x direction for all 
values of w, corresponding to the fact that the solid black 
curve in Figs. 4(a), 4(b), 4(c) falls inside the shaded re- 
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(a) E = 50, uj = 1 



(b) E = 17, u = 1 
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(c) E = 50, u> = 
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(d) E = 17, uj = 
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(e) E = 50, uj = -1 



(f) E =17,ui 



Figure 5: Boundaries of motion for different values of the 
energy and the u parameter, for string starting at rest at 
r = th + 1.25 and 6 = tt/2 + 0.02. The grey area represent 
the region inside the event horizon, the dashed black curve 
the boundary of the ergoregion and the dashed red line the 
boundaries of motion. The form of the boundaries is exactly 
what one could have predicted by using Figs. 4(a)-4(c). 



gion when E = 50. Additionally, for uj — 1,0 there is no 
boundary along the y direction, whereas for uj = — 1 there 
is such a boundary, corresponding to the fact that in the 
latter case the E — 50 point on the solid black curve lies 
below the dashed red line in Fig. 4(c). In Fig. 6, one can 
see the region closer to the black hole: for the uj = 0, — 1 
cases, there exists a third boundary, which is an outer 
boundary. This is in accordance with the form of Ef, in 
Figs. 3(b) and 3(c) and highlights the subtlety discussed 
above right after outlining the various regions of Fig. 4. 
For E — 17 on the other hand, only the uj = 1 choice al- 
lows motion which is unbounded in the y direction. For 
uj = the string is trapped in a toroidal volume, and for 
u> = — 1 there is only an outer boundary. Fig. 7 is similar 
to Fig. 5, but for r — nsco- (We do not plot the u = — 1 



(a) E = 50, uj = 




(b) E = 50, uj ■■ 



Figure 6: Magnified segments of Fig. 5(c) and 5(e) which 
demonstrate the existence of yet a third boundary of motion 
for uj = 0, — 1. This acts as a second outer boundary, therefore 
forbidding any string starting from inside the region bounded 
by the innermost red dotted curve to ever leave this region. 



case since there is no energy for which there will be an 
internal boundary for the motion.) In this case, the inner 
boundary of motion lies within the ergoregion. 

Before closing this section it is worth stressing once 
more the similarity between the contour plots of the 
boundaries of motion for a Newtonian elastic spinning 
ring presented in Fig. 1 and the boundaries of motion for 
the relativistic strings presented in this section. 



3. Trajectories 

In this subsection we look at examples of individ- 
ual trajectories, and in particular we exhibit processes 
whereby a string is ejected along the axis, with some 
of the initial internal energy of the string converted to 
translational kinetic energy. To this end, we employed 
numerical integration with Mathematica 6.0 to solve the 
equations of motion. As already mentioned we consider 
here two different initial conditions: one for which the 
string starts outside the ergoregion, and one for which 
the string starts at r — nscOi which lies well within 
the ergoregion (points P\ and P2 in Fig. 2 respectively). 
Starting at rest from these initial positions, there are 
three types of trajectories, as indicted by the types of 
boundaries of motion: those that are trapped inside an 
outer boundary with no inner boundary and simply fall 
in across the horizon, those that are trapped in some 
toroidal volume, and those that can reach spatial infin- 
ity along the direction of the axis. The trajectories of 
the first type display no interesting features so we do not 
present any explicit examples of them here. 7 



7 Note that one can have trajectories of this type even if there 
exists an inner boundary, provided that an outer boundary at 
smaller x's exists as well, as discussed earlier. This is the case, 
for instance, when the boundaries of motion are those of Fig. 6(a) 
or Fig. 5(d) and the string start at rest at r = + 1.25. 
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(a) E = 38, w = 1 



(b) E = 8, u> = 1 
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(c) £ = 38, w = 1 



(d) B = 8, w = 1 
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(e) £ = 38, w = 



(f) E = 8, u = 



Figure 7: Boundaries of motion for different values of the 
energy and the u parameter, for string starting at rest at 
r = nsco and 8 = tt/2 + 0.02. The form of the boundaries 
is exactly what one could have predicted by using Figs. 4(d)- 
4(f). Figs. 7(c) and 7(d) focus on the region close to the 
black hole of Figs. 7(a) and 7(b) respectively. Clearly, the 
inner boundary of motion lies within the ergoregion. 



We begin by discussing bound trajectories. The en- 
ergies we considered previously do not lead to confined 
motion for a string starting at P\. This can be seen in 
Fig. 5: For the cases where there can be confined mo- 
tion, namely Figs. 5(d) and 5(e), our initial position lies 
on the innermost outer boundary. Thus the string would 
fall straight into the black hole. However, for other ener- 
gies there can be confined motion for a string starting at 
Pi. Also, according to Fig. 7(b) we have already found 
an example where the string's motion is confined when 
it starts from P2, i.e. inside the ergoregion. In Fig. 8 
we present the motion of a string starting from Pi with 
E = 9 and that of a string starting from P 2 with E = 8 
for qualitative comparison. Both cases refer to u> = 1 





(a) E = 9, r = r h + 1.25 



(b) E = 9, r = r h + 1.25 




10 15 20 25 1 




(c) E = 8, r = r lsco 



(d) E ■■ 



r = nsco 



Figure 8: Confined trajectories of the string for different ini- 
tial positions and different energies. The blue continuous line 
corresponds to the trajectory of a point on the string ((f> di- 
mension suppressed). No qualitative difference seems to exist 
between the motion of a string starting inside the ergoregion 
and one starting outside. In both case the motion seems to 
follow no clear pattern. 



(co-rotation). In both cases the strings move somewhat 
randomly, without following any clear pattern. Fig. 8(b) 
and 8(d) focus on the area inside and around the ergore- 
gion. The string is sometimes pulled by gravity rather 
sharply back toward the equatorial plane after reflecting 
from the inner boundary. We have monitored the ac- 
curacy by running the code backward in time from the 
point the forward run ended, and checking how closely 
this backward run approaches the initial conditions. In 
the plots presented we stop the simulation before the ac- 
curacy becomes lower than 10 -4 . 8 

Before proceeding further to trajectories that can reach 
spatial infinity, there is another interesting feature of con- 
fined motion which should not be missed. In Fig. 9 we 
present confined trajectories starting from Pi with E = 8 



8 Using only the default settings of Mathematica, we were able 
to run the simulation for r ~ 150 — 200 for trajectories start- 
ing from Pi and for r ~ 50 — 100 for trajectories starting from 
P2 (provided that the string did not end up in the inside the 
horizon before that time) before dropping below the required ac- 
curacy. Remarkably, for the "focused" trajectories, such as those 
in Fig. 9, the accuracy was still two orders of magnitude higher 
than our minimum requirement at r = 500. Note that r is the 
coordinate time on the worldshcct. 
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(a) E = 8, r = r h + 1.25 (b) E = 8, r = r h + 1.25 





(c) E = 5.5, r = nsco 



(d) E = 5.5, r = nsco 



Figure 9: "Focusing" of the trajectories. In this case a pattern 
seems to be arising. In Figs. 9(b) and 9(d) we have distorted 
the aspect ratio in order to make this manifest. 



and from P2 with E — 5.5, always for u) = 1. There is 
evidently some remarkable "focusing" of the motion close 
to the equatorial plane in both cases, and a very regu- 
lar pattern is generated. In Figs. 9(b) and 9(d) we have 
distorted the aspect ratio in order to make manifest the 
pattern that arises in the motion. We do not under- 
stand why these cases are so different from those shown 
in Fig. 8. The focussing cannot be explained by consider- 
ing only the boundary of allowed motion (dotted curve) 
as inferred using the effective potential. The appearance 
of a Lissajous-like pattern indicates that when the string 
stays close to the equatorial plane the motion becomes 
similar to complex periodic oscillation. 

Finally, we discuss trajectories that can reach spatial 
infinity. 9 In Fig. 10(a) we have plotted the trajectory 
of a string starting from P\ with E = 17 and u> = 1. 
After a couple of oscillations around the equatorial plane, 
which arc exhibited more clearly in Fig. 10(b), the string 
escapes along the z direction. The behavior of a string 
starting from Pi with E — 38 and u) = 1 is quite similar. 
In both cases, there is some conversion of internal energy 
of the string into kinetic energy along the y axis. 

Another interesting class of non-confined trajectories 
are the scattering trajectories found by Larsen in Ref. [1]. 
An example of such a trajectory for a string with E = 




(a) E = 17, r = r h + 1.25 (b) E = 17, r = r h + 1.25 




10 20 30 40 50 



(c) E = 38, r = r ISCO 



(d) E = 38, r = r ISCO 



9 These trajectories are similar to the time-reversed version of 
Larsen's adiabatic capture trajectories found in Ref. [1]. 



Figure 10: Escape trajectories for strings of different energy 
and angular momentum, start at Pi and P2 respectively. Af- 
ter some oscilations around the equatorial plane the strings 
escape to infinity along the axis of symmetry. Clearly there 
is some conversion of internal to kinetic energy. 



11.2929, u> = 1, and initial conditions r = 200, #0 = 
3.11, r = -3.4 and O = -0.007 is shown in Fig. 11(a). 
The string comes in along the axis and after a couple of 
bounces around the black hole, which can be better ap- 
preciated in Fig. 11(b), escapes to infinity out the other 
side. The scattering in this case reduces the mean radius 
in the x direction, signaling a conversion of internal en- 
ergy to kinetic energy along the y direction, as explained 
below in the next subsection, VB4. 

In Figs. 11(c) and 11(d) we show a case where the 
string actually falls into the black hole instead if getting 
scattered. For this case, E = 11.8417, ui — 1, and the 
initial conditions are ro = 200, #0 = 3.11, ro = —2.4 and 
6q = —0.007. In this case the string has more energy than 
a string at rest on the horizon and, hence, its motion has 
no inner boundary close to the equatorial plane (cf. the 
discussion in section VB2 about Figs. 3(a) and 4). 



4- Ejection energy 

As noted, on some trajectories internal energy of the 
string can be converted to translational kinetic energy 
along the axis. To characterize this energy we may as- 
sociate a relativistic boost factor V = E/Eq with the 
asymptotic string, where Eq is the energy of the string in 
its center of mass frame at infinity. In the center of mass 
frame, all velocities vanish at the inner and outer bound- 



12 





(a) E = 11.3, scattering (b) E = 11.3, scattering 





10 15 20 



(c) E = 11.8, collapse 



(d) E = 11.8, collapse 



Figure 11: Scattering and collapse on the black hole of strings 
coming from spatial infinity. The scattering can lead to energy 
conversion, as can be seen from Fig. 11(a) by the fact that 
the average radius of oscillation along the x direction has been 
slightly reduced after the scattering. On the other hand, if 
the inner boundary of motion lies within the horizon close to 
the equatorial plane, then instead of a scattering there can be 
a collapse, see Figs. 11(c) and 11(d). 



can be produced. This statement is based both on exam- 
ining cases (not explicitly reported here) for which a = 0, 
as well as by comparing the behaviors for different values 
of uj. The dynamics only depends on uj when a is non- 
zero, but the effect of changing uj can also generally be 
reproduced by changing other quantities. 

The one exception to this statement concerns the issue 
of negative Killing energy states. In the presence of black 
hole spin, there is an ergoregion, inside of which such 
states exist for point particles. The same is true for the 
string, though the mechanism is a bit different. 

For a particle, negative angular momentum orbits in 
the ergoregion have negative Killing energy. Longitudi- 
nal motion of a string is meaningless, so it is not the "or- 
bital" motion, but rather the current on the worldsheet 
that can access negative energy states. In particular, the 
second term on the right hand side of the energy (36), 
i.e. —gt ( j 1 Ti' JT — \g t j > /g c / >c /,\L/2ir, is negative there when 
the angular momentum L = 4itujJ 2 /(l + uj 2 ) is negative. 
The first term in eq. (36) vanishes at the horizon, so the 
energy is certainly negative there when L < 0. The neg- 
ative energy states extend out a certain distance (which 
depends on the values of uj and J) from the horizon into 
the ergoregion. For example, one can see in Fig. 4(f) the 
negative energy states at nsco in the case uj = — 1, ex- 
cept for the smallest values of J. (For the limiting case 
of zero tension and uj = — 1, there are negative energy 
states everywhere in the ergoregion.) 



VI. RELEVANCE TO ASTROPHYSICS? 



aries p — pi >0 of oscillation, so we may use (43) to obtain 
E = p+J 2 /p. Eliminating J we find E = pi+p . Thus, 
using the present notation x for the cylindrical radius, we 
have 

r = — E — r . (so) 

X% T Xq 

With this relation the T factors corresponding to the pro- 
cesses shown in the previous examples can be computed. 
For example, for the E — 17 case in Fig. 10 we find 
r w 1.01, or v/c » 0.15, and for the E — 38 case we find 
r w 1.09, or v/ck, 0.39. We expect that these mildly rel- 
ativistic boost factors are typical of what one can expect 
in generic cases. 



5. The role of black hole spin 

The black hole spin enters the string dynamics via the 
a dependence of the various metric coefficients. While 
the value of a, certainly affects the boundaries of motion 
and dynamics when other quantities are held fixed, we 
do not see any systematic qualitative effect that can be 
attributed uniquely to the spin. That is, by varying the 
current, energy, and location of the string similar effects 



As stated in the Introduction, although our primary 
aim in this paper was to study the motion of a current 
carrying strings in a Kerr background as an interesting 
dynamical system, our investigation was also loosely mo- 
tivated by the problem of the production of collimated 
jets in astrophysics. We should therefore like to com- 
ment on what might be the relevance of our findings to 
astrophysical systems. 

First of all, the current carried by the string in the 
model need not represent an electric current, but it can 
just serve to model the role of angular momentum. In 
this sense, a current carrying string can serve as an ide- 
alization of some MHD matter configuration with tension 
and angular momentum, as mentioned in the Introduc- 
tion. Simplistic as it may be, this allows for the identifica- 
tion of a mechanism that might conceivably be relevant: 
the tension sets a barrier for the motion that keeps the 
string within a certain radius from the axis of symmetry, 
while the angular momentum sets a barrier that keeps 
the string outside a certain radius. Gravity deforms these 
barriers creating either a closed region in which the string 
is trapped, or a "restricted path" that leads to motion 
along the symmetry axis. If some physical process were 
to "feed" such a system with strings close to the equato- 
rial plane, then under the right circumstances one might 
end up with a very well collimated stream of strings mov- 
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ing out along the axis of symmetry. Of course, the ax- 
isymmetry we are imposing might be too restrictive, and 
in particular, even if a string started out approximately 
axisymmetric, it might wobble or twist and reconnect. 
Perhaps the tension could act to suppress such effects, 
but this requires further investigation. 

Still, one may speculate that the general flavor of the 
mechanism exhibited by the strings might generically be 
present in a class of astrophysical systems. For instance, 
given a system with an accretion disk, plasma leaving 
the inner edge of the disk would carry some angular mo- 
mentum, and plasma flux tubes do possess a significant 
tension (which however is not constant but grows in pro- 
portion to the length) . It might even be that differential 
rotation can stretch such a tube, storing in it a large 
amount of electromagnetic energy before launching it to- 
wards the axis. External magnetic fields, which have 
been left out of consideration in our initial study, might 
also play an important role in the subsequent dynam- 
ics. It remains to be seen whether aspects of this picture 
could actually be astrophysically relevant. 

Regarding the energetics, it should be noted that the 
mechanism described here does not lead to significant 
(ultra-relativistic) acceleration of the strings. While in- 
ternal energy of the string can be converted to kinetic 
energy of motion along the axis of symmetry, the result- 
ing T factors in the examples we studied are rather close 
to unity, the largest velocity we reported being 0.39c. 
Therefore, even if the mechanism we have found is some- 
how relevant to collimation, it does not appear able to 
explain the acceleration of relativistic jets. The acceler- 
ation would have to occur due to the presence of elec- 
tromagnetic fields, or to some other mechanism. In par- 
ticular, the process discussed here does not tap into the 
rotational energy of the black hole, which is a prime can- 
didate for the source of ultra-relativistic jet energy. 

On the contrary, the collimation effect we find persists 
even at the Newtonian level, which may be advantageous: 
collimated jets seem to be present not only around black 
holes, but also around objects with significantly weaker 
gravitational fields, such as accreting young stars for in- 
stance [3]. Therefore, if there is a common mechanism 
for the collimation, it should not depend on relativistic 
effects. Hence, the effects described here could perhaps 



be part of a more involved process that collimates and 
drives astrophysical jets. 



VII. CONCLUSIONS 

We have studied the motion of relativistic, current car- 
rying strings moving axisymmetrically on the background 
of a Kerr black hole. First the equations of motion and 
conserved quantities in the conformal gauge were found. 
Next, we determined how the energy of a string that is 
instantaneously at rest depends on the current and the 
position of the string. Contours of this quantity deter- 
mine boundaries of the possible motion of a string with 
that energy and current. We compared this analysis with 
that for a Newtonian elastic ring moving around a point 
mass, and found that system to be qualitatively similar. 

By considering a number of examples and various plots 
of the relevant functions, the possible types of motion 
were mapped out. Regions of parameters for which the 
string falls into the black hole, or is trapped in a toroidal 
volume, or can escape to infinity, were identified. After 
this general analysis, we examined representative trajec- 
tories found by numerical integration, illustrating vari- 
ous interesting behaviors. In particular, we found that 
a string can start out at rest near the equatorial plane 
and, after bouncing around, be ejected out along the axis, 
some of its internal (elastic or rotational kinetic) energy 
having been transformed into translational kinetic en- 
ergy The resulting velocity can be an order unity frac- 
tion of the speed of light. 

Finally, we addressed the question of possible astro- 
physical significance of this system as a simple model of 
MHD plasma flux tubes, which might conceivably play a 
role in jet formation and collimation. 



Acknowledgments 

We thank Vijay Kaul for collaboration in the early 
stages of this research. This work was supported by 
the National Science Foundation under grant PHYS- 
0601800. 



[1] A. L. Larsen, "Chaotic string capture by black hole," 
Class. Quant. Grav. 11, 1201 (1994) [arXiv:hep- 
th/9309086]. 

[2] A. V. Frolov and A. L. Larsen, "Chaotic scattering and 
capture of strings by a black hole," Class. Quant. Grav. 
16, 3717 (1999) [arXiv:gr-qc/9908039]. 

[3] M. Livio, "Astrophysical jets: a phenomenological exam- 
ination of acceleration and collimation," Phys. Rep. 311, 
225 (1999). 

[4] M. Christensson and M. Hindmarsh, "Magnetic fields in 
the early universe in the string approach to MHD," Phys. 



Rev. D 60, 063001 (1999) [arXiv:astro-ph/9904358]. 
[5] V. S. Semenov, S. A. Dyadechkin, I. B. Ivanov and 

H. K. Biernat, "Energy confinement for a relativistic 

magnetic flux tube in the ergosphere of a Kerr black 

hole," Phys. Scripta 65, 13 (2002). 
[6] V. Semenov, S. Dyadechkin and B. Punsly, Science 305, 

978 (2004) [arXiv:astro-ph/0408371]. 
[7] H. C. Spruit, "Equations for thin flux tubes in ideal 

MHD," Astron. Astrophys. 102, 129 (1981). 
[8] A. L. Larsen, "Dynamics of cosmic strings and springs: 

A Covariant formulation," Class. Quant. Grav. 10, 1541 



14 



(1993) [arXiv:hep-th/9304086]. 
[9] A. Villcnkin and E. P. S. Shellard, Cosmic Strings and 
Other Topological Defects, (Cambridge University Press, 
Cambridge, 1994). 
[10] M. B. Hindmarsh and T. W. B. Kibble, "Cosmic 
strings," Rept. Prog. Phys. 58, 477 (1995) [arXiv:hep- 
ph/9411342]. 

[11] B. Carter, P. Peter and A. Gangui, "Avoidance of col- 
lapse by circular current-carrying cosmic string loops," 
Phys. Rev. D 55, 4647 (1997) [arXiv:hep-ph/9609401]. 

[12] B. Carter, "Brane dynamics for treatment of cosmic 
strings and vortons," arXiv:hep-th/9705172. 

[13] A. L. Larsen, "Charged circular test string in external 
magnetic fields," Mod. Phys. Lett. A 7, 2913 (1992). 

[14] A. L. Larsen, "Charge current carrying circular strings in 
arbitrary backgrounds," Phys. Lett. A 170, 174 (1992). 

[15] S. A. Dyadechkin, V. S. Semenov, H. K. Biernat and 
T. Penz, "Comparison Of Magnetic Flux Tube And Cos- 
mic String Behavior In Kerr Metric," Adv. Space Res. 



42, 565 (2008). 

[16] F. Dubath, M. Sakellariadou and C. M. Viallet, "Scat- 
tering of cosmic strings by black holes: Loop forma- 
tion," Int. J. Mod. Phys. D 16, 1311 (2007) [arXiv:gr- 
qc/0609089]. 

[17] M. Snajdr and V. P. Frolov, "Capture and critical 
scattering of a long cosmic string by a rotating black 
hole," Class. Quant. Grav. 20, 1303 (2003) [arXiv:gr- 
qc/0211018]. 

[18] M. Snajdr, V. P. Frolov and J. P. DeVilliers, "Scattering 
of a long cosmic string by a rotating black hole," Class. 
Quant. Grav. 19, 5987 (2002) [arXiv:gr-qc/0208009]. 

[19] E. Witten, "Superconducting Strings," Nucl. Phys. B 
249, 557 (1985). 

[20] J. M. Bardeen, W. H. Press and S. A. Teukolsky, "Ro- 
tating Black Holes: Locally Nonrotating Frames, Energy 
Extraction, And Scalar Synchrotron Radiation," Astro- 
phys. J. 178, 347 (1972). 



